D D C 


♦Research  supported  by  NRC  Grants  lA-3973,  A-3990,  and  A-9239  at  the 
University  of  British  Columbia,  and  by  U.S.  Army  research  Contract 
#DAAG-29-76-G-0210. 

♦♦Department  of  Mathematics,  Univeralty  of  Utah,  Salt  Lake  City,  Utah  84112. 


A . • %.  m -4k 


A SINC-GALERKIN  METHOD  OF  SOLUTION  OF  BOUNDARY  VALUE  PROBLEMS 

Frank  Stenger,  Department  of  Mathematics, 
University  of  British  Columbia 
Vancouver,  B.C. 

Abstract 


This  paper  illustrates  the  application  of  a "Sinc-Galerkin"  method 
to  the  approximate  solution  of  linear  and  nonlinear  second  order 
ordinary  differential  equations,  and  to  the  approximate  solution  of 
some  linear  elliptic  and  parabolic  partial  differential  equations  in 


the  plane.  The  method  is  based  on  approximating  functions  and  their 

derivatives  by  use  of  the  Whittaker  cardinal  function.  The  DE  is 

reduced  to  a system  of  algebraic  equations  via  new  accurate  explicit 

approximations  of  the  inner  products,  the  evaluation  of  which  does  not 

require  any  numerical  integration.  Using  n function  evaluations 

the  error  in  the  final  approximation  to  the  solution  of  the  DE  is 
l/2d 

_cn 

0(e  ) where  c is  Independent  of  n,  and  d denotes  the 

dimension  of  the  region  on  which  the  DE  is  defined.  This  rate  of 


convergence  is  optimal  in  the  class  of  n-polnt  methods  which  assume 
that  the  solution  is  analytic  in  the  interior  of  the  Interval,  and 
which  ignore  possible  singularities  of  the  solution  at  the  end-points 


of  the  Interval. 
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1.  Introduction  and  Summary 

The  function  slnc(x)  Is  defined  on  the  real  line  by 
(1.1)  sinc(x)  - 

The  Whittaker  cardinal  function  of  an  arbitrary  function  f Is  defined 
for  any  h > 0 by 


(1.2)  C(f ,h,x)  - l f(kh)sinc[£pl],  h > 0, 

k— <» 


whenever  this  series  converges. 

The  approximation  of  f using  a finite  number  of  terms  of  (1.2) 

has  been  extensively  studied.  The  paper  [8]  contains  a review  of  the 

properties  of  C(f,h,x)  which  were  discovered  by  E.T.  Whittaker  [15], 

J.M.  Whittaker  [16],  Hartly  [5],  Nyquist  [9]  and  Shannon  [12].  In  [13] 

new  approximations  are  derived  by  means  of  C(f,h,x),  for  interpolating, 

integrating  and  approximating  the  Fourier  (over  (-«,<•)  only)  and  Hilbert 

transforms  over  (-“,“) , (0,«)  and  (-1,1)  . In  [fl  the  function 

C(f,h,x)  is  used  to  obtain  formulas  for  approximating  the  derivatives 

of  functions  over  (-*•,"),  (0,»)  and  (-1,1)  . 

In  the  present  paper  we  use  results  of  [7,13]  to  derive  basis 

functions  for  Galerkln  schemes  of  solving  second  order  problems, 

and  we  derive  explicit  and  highly  accurate  expressions  for  inner 

d2  d 

products  such  as  (f  , *k),  (f  ^ , ♦fc),  (fu,  *k)  . All  of 

these  are  expressed  in  terms  of  the  function  values  of  u,  and  not 
tha  derivatives  of  u . We  then  study  the  application  of  the  derived 


f 


[ 

I 
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approx last ions  on  Che  approximate  solution  of  some  ordinary  and  partial 
differential  equatlona,  via  the  Galerkln  method.  The  combined  method 
thus  ylelde  a system  of  algebraic  equations,  without  the  uae  of  any 
quadrature. 

Let  us  briefly  compare  the  present  method  of  approximate  solution 
of  linear  differential  equations  with  currently  popular  finite  difference 
methods,  or  with  finite  element  methods  that  use  piecewise  linear 
elements.  The  finite  difference  or  finite  element  methods  lead  to  a 
sparse  system  of  equations.  The  use  of  n solution  evaluations 
usually  leads  to  a linear  system  of  n algebraic  equations  having 
non-singular  coefficient  matrix.  The  error  in  the  resultlon  approxi- 

'utlon  is  0(n  P) , where  p is  usually  1 or  2 . The  "Sinc- 
Galerkin"  method  of  the  present  paper  also  leads  to  a system  of  order 
n on  the  basis  of  n solution  evaluations.  This  system  has  a 
nonsingular  full  matrix.  The  error  in  the 

-cn1/(2d) 

resulting  approximate  solution  is  0(e  ),  where  d denotes 

the  dimension  of  the  problem.  The  advantage  of  the  present  method  is 
that  due  to  its  rapid  convergence  it  does  not  require  the  solution  of 
a very  large  system  of  equations  in  order  to  achieve  a dejiyed  accuracy, 
if  more  than  two  significant  figures  of  accuracy  are  required  in  the 
approximate  solution.  In  addition,  the  rate  of  convergence  of  the 
preaent  method  is  the  aame,  regardleaa  of  posalble  singularities  of 
the  solution  of  an  equation  on  the  boundary  of  the  region. 


The  approximate  methods  of  [7,13]  have  previously  been  affectively 
applied  to  the  approximate  aolutlon  of  integral  equations  via  Galerkln- 


type  methods  In  [2,10,11].  In  [6]  an  effective  Galerkin-type  method 
Is  derived  which  uses  approximations  derived  in  [13]  to  obtain  an 
approximate  solution  to  the  problem 

(1.3)  y"  - y - y3/-2  , y(0)  - y(-)  - o 

via  the  minimization  of  a certain  nonlinear  functional.  In  all  of  the 

-cn* 

above  cases  the  error  of  an  n-point  approximate  solution  is  0(e  ) . 

In  Sec. 2 of  the  present  paper  we  review  the  relevant  known 

approximation  properties  of  Whittaker's  cardinal  function,  and  we  then 

use  these  to  derive  explicit  approximate  inner  products,  in  general  as 

well  as  for  the  important  special  cases  of  the  intervals  [0,1],  [-1,1], 

[0,«»]  and  [-*,■*]  . In  Sec. 3 we  illustrate  the  application  of  the 

previously  derived  formulas  to  the  approximate  solution  of  some  simple 
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"model"  problems,  such  as  u"  ■ -2,  u"  ■ u - u /x  , ut  ■ and 

u + u ■ f,  with  appropriate  boundary  conditions.  In  Sec. 4 we 

XX  yy  -cn1/(2d) 

carry  out  an  error  analysis,  proving  the  0(e  ) rate  of 

convergence  referred  to  above. 

It  is  worthwhile,  for  purposes  of  better  understanding  of  the 
results  of  the  paper,  to  consider  the  solution  of  the  one-dimensional 
problems 

L(f ) (x)  - f"(x)  + q(f (x) , f'(x).x)  - 0,  0 < x < 1 

(1.4) 

f (0)  - f (1)  - 0 

or 


L(f)(x)  - rn(x)f"(x)  + r. (x)f ' (x)  + r,(x)f(x)  • r_(x),  0 < x < 1 
(1.5)  u 1 J 

f(0)  - f(l)  - 0 





We  seek  an  approximate  solution  of  (1.4)  or  (1.5)  of  the  form 


N 

(1.6)  f (x)  ■ J f.  S(k  hW(x) 

k— N 


where 


sin[£(x-kh)] 

(1.7)  S(k,h)(x)  - j2 . *(x)  - log  TIT  • 

J(x-kh) 

We  introduce  a convenient  function  g which  plays  the  role  of  a weight 
function,  set 

(1.8)  SR(x)  - g(x)S(k,h)«$(x) 

and  reduce  the  solution  of  (1.4)  or  (1.5)  to  the  solution  of  the  Galerkin 
system 

1 

(1.9)  (L(f),Sk)  - j L(f)(x)Sk(x)dx  - 0,  k - - N,  - N+l N . 

The  results  of  thm.  2.11  together  with  Eq.  (2.54)  may  then  be  used  to 
get  an  explicit  approximation  of  (1.9)  for  the  case  of  (1.5),  whereas 
the  combination  of  Thm.  2.7,  2.11  and  Eqs.  (2.54)  may  be  used  to  get 
an  explicit  approximation  of  (1.9)  for  the  case  of  (1.4). 

The  author  is  grateful  to  A.  Adler  and  J.  Varah  at  the  U.B.C. 
for  valuable  discussions  lvolvlng  Thm.  2.12  of  this  paper.  The  author 
also  wishes  to  thank  Bob  Burke  of  the  Unlv.  of  Utah  for  carrying  out  the 
computations  in  Sec.  3. 


2 . Preliminaries  and  Fundamentals 


In  Chls  section  we  shall  recall  some  known  properties  [8]  and 
derive  some  new  properties  of  Whittaker's  cardinal  function,  which  we 
shall  require  in  this  paper. 

Definition  2.1.  Let  R denote  the  real  line,  C the  complex  plane,  and 
let  B(h)  denote  the  family  of  all  functions  defined  on  C that  are 
entire,  such  that  f (R)  and  such  that 

(2.1)  |f(z)|  < Ce1,!yl/h,  z-x  + iy^C, 
for  some  constant  C . Set 

(2.2)  S(j,h)(x)  - sincL2^] 
and 

(2.3) 


6(n)  _ s(n)(j,l)(k)  - <^)nS(j,D(x)|x.k 


In  particular,  we  have 


(2.4) 


(0)  . I1  u 3 

- k 

Jk 

^0  if  j 

i k 

;(1)  ■ 

. f i; 1 

- k 

Jk 

H-l)  j 

if 

k-J 

,(2) 

Jk 


-i 


- 5 if  J " k 


-2(-Dk~j 

(k-j)2 


if  j i k , 


Theorem  2.2  [8]:  Let  f £ B(h)  . Then  f(z)  - C(f,h,z)  . Moreover 

ir_ 

(2.5)  f(z)  - | g(t)elztdt  for  some  g f.  L2(-  £ , £ ) ; 


(2.6) 


(t)dt  ; 


f °° 

(2.7)  |f(x)|2dx  -hi  | f (kh) | 

J R k <■-<*> 


and  the  sequence  (h"*S(k,h) is  therefore  a complete  orthonormal 
sequence  In  B(h); 


(2.8) 


f 6 B(h)a^f'  & B(h) 


Theorem  2.3:  Let  6:11/  be  defined  as  in  (2.3).  Then 


(2.9)  | (sinc[x~^h]  )^n^sinc[xh— 1 dx  - h1  n 6^  , 


n — 0|1|2|...  . 


Proof:  Let  us  set 


(2.10) 


f(t)  - S(j,h)(t) 


and  let  us  note  that  f G B(h)  . By  Eq.  (2.8)  it  thus  follows  that 
f(n)  6 B(h) , n - 0,1,2,...  . Eq.  (2.9)  thus  follows  by  taking 
f - S(J,h)(n)  in  (2.6),  and  noting  by  (2.3)  that 


(2.11) 


S(n)(j,h)(kh)  - h‘n  6™  . 


Definition  2. A:  Let  d > 0,  and  let  B(Pp  denote  the  family  of  all 
functions  f that  are  analytic  in 

(2.12)  V'd  - (z  - x + ly  : |y|  < d), 
such  that 

rd 

(2.13)  | f (x+iy) ) dy  -*•  0 as  x -*•  ±® 

-d 

and  such  that  N(f,P^)  < ®,  where 

(2.14)  N(f,PI)  - lira  {I  |f(x+iy)|dx  + f |f(x-iy)|dx} 

y-*d~  Jr  Jr 

Theorem  2.5  [13]:  Let  h and  d be  positive,  let  f £ B(Pp,  and 
let  c(f)  be  defined  by 

(2.15)  e(f ) (x)  ■ f(x)  - C(f,h,x),  x £ R . 

Then 

/->  i*\  ,.w  \ 9ln(mc/h)  f , f (t-id~) f (t+id~) ,. 

(2.16)  e(f)(x)  - 2iri  j ^ (t-x-id)sin[ t-id)ir/h)  “ (t-x+id)sinlt+id)ir/h] 

Moreover 

N(f.P') 

(2.17)  ^P'E(f)(x)'  - 2nd  sinh(itd7h)  * 

Definition  2.6:  Let  P be  a simply-connected  domain  in  the  complex 

plane  C,  and  let  PJ  be  defined  as  in  (2.12).  Let  ♦ be  i conformal 

map  of  P onto  PI,  and  let  V * ♦*1  denote  the  inverse  map.  Let 
d 

• m ii(-m)  and  b - 4>(m)  j*  a be  boundary  points  of  P,  and  let  us  take 


(2.18) 


r « (v  £ P : w - ip(x)  , -•  i x <_  •»}  . 


Let  B(P)  denote  the  family  of  all  functions  that  are  analytic  In  P, 


such  that  for  u real 


(2.19) 


where 


| f (z)dz  | 0 as  u -*■  ±°» 

JtKL+u) 


(2.20) 


L*  {ly  : -d  < y < d}, 


and  such  that 


(2.21)  N(f,P)  - 11m  Inf  f |f(z)dz|  < •» 

OSP.CCP  'C 


(Note  that  If  f 6B(P),  then  f • ip  tB^).)  Set 


(2.22) 


- ip(kh)  , k - 0,11,12, . . . , 


and  let  g be  a function  which  is  analytic  in  P,  which  plays  the  role  of  a 
weight  function  In  the  inner  products,  and  whose  properties  we  shall  determine 
In  the  sequel.  Finally,  we  set 


(2.23) 


Sj(z)  « g(z)sinc(^~^]  - g(z)S(J  ,h)o$(z)  . 


The  following  result  was  established  in  [7], 


Theorem  2.7;  Let  m be  a nonnegative  Integer,  and  let  f^'/g  £ B(P)  . 
Let  there  exist  positive  constants  a , C©  depending  only  on  m,  d 
and  g,  depending  only  on  m and  g,  and  Cj  depending  only  on 

m,  g and  f,  such  that 


(2.24) 


< C,  e"°I*(x)|  for  all  x 6 T 
g(x)  - 2 


(2.25) 


(2-)nS,  (x)  < C,  h for  all  x £ T 

dx  k “1 


(2  26)  I (j_)n(B(x)a1n[x^)/hljl 

(2-26)  |ldx)  \ ♦(z)  - ♦(*)  J|-  0 


n-0,1, . . . ,m. 


h for  all  x € T,  z £ 3DJ 


Then  there  exists  a constant  K depending  only  on  m,d,a,g  and  f such 

jtf 

that  if  h - [ttd/(aN)P  then 


N f(x.) 


(2.27) 


|f<n,(x)  - l ;(x)|  < K N exp[-(*daNK] 

j— N gUj'  J 


for  all  x £ f,  and  for  n ■ 0,1,. . . ,m 


Theorem  2.8  [13]:  If  f fB(P),  then  the  identity 


« f(x.) 


(2.28) 


TTt\  - l TT7^~\  S(j,h)H 
*'<*>  jL  ^j) 


sin[x^(x) 

2iri 


f f(z)  dz  

j [$(z)-<t>(x)]sin[ir<t>(z)/h] 


is  valid  for  all  x fF  . Moreover 


r - f(x.) 

I ' h jL 

r J J 

^ r exp[**-^Z-  sgn  Im  $(z)] 

2L  sin[i*(z)] 


(2.29) 


f (z)dz 


in[-  $(z)  1 


The  results  of  this  theorem  may  be  conveniently  combined  with  those  of 

the  formulas  obtained  above,  to  yield  explicit  approximate  expressions 

for  inner  products.  The  results  of  the  following  lemma  are  useful  for 

*Here  and  henceforth  f(z)dz  is  defined  by  lim  inf  f(z)dz,  for  any 

L O30,CCt>Jc 

f € B(P)  . 3P 


bounding  the  error  of  these  approximate  expressions. 

Lemma  2.9:  If  |lm  z|  - d > 0 and  If  k is  an  integer,  then 


(2.30) 


sincf (z-kh) /h] 
sin(nz/h) 


i Cl(h’d>  5 2 Sd 


(2.31) 


(d/dz){slnc[ (z-kh) /h] } 
sin(irz/h) 


< C (h,d)  3 d-Kh/n)tanh(xd/h)  . 

2d  tanh(rtd/h) 


(2.32) 


dz 


2<sinc[ (z-kh)/h] } 


sin(x/h) 


< C (h  d)  = t(2h/x)+TT  d/h]d  tanh(*d/h)+2d 
3 2 d3tanh(ird/h) 


(2.33) 


|sin(irz/h)|  ^ sinh(ird/h) , |cos(irz/h)|  4 cosh(ird/h)  • 


Proof : We  shall  only  prove  (2.31),  since  the  proofs  of  the  remaining 
cases  are  similar,  and  we  omit  them.  We  have 


— { s(k  h)(z)}  - cos[x(z-kh)/h]  _ h sin[x(z-kh)/h] 

dz  * z-kh  n , ...2 

(z-kh) 


Now  if  | Imz | ■ d,  then  | r— kh | ^ d,  | cos [ir( z-kh) ] | <,  cosh(ird/h)  and 
|sin(irz/h)|  ^ sinh(itd/h)  ; hence 


w 

siri(xz/h) 


“ d tanh(ird/h)  + 


C2(h,d)  . 


Theorem  2.10:  Let  6^"^  be  defined  as  in  (2.3)  and  (2.4),  let  C,(h,d) 
be  defined  as  in  Lemma  2.9,  let  x^  be  defined  as  in  (2.22),  as  in 

(2.23),  set  F^  - F(x^)  for  an  arbitrary  function  F,  and  1st  r and 
f be  functions  which  are  analytic  in  V . 


(a)  Let  rfg  6 B(P)  . Then 


(2.34) 


fkrk8k 


r(x)f(x)S  (x)dx  - h *.*■  - < C (h,d)N(rfg,P)e 

r 9k 


-xd/h 


(b)  Let  trfg/$](x)  ■*>  0 as  x a and  aa  x -*■  b along  T,  and  let 
(rg)'f  and  rg$'f  ^B(P)  . Then 


(2.35)  r(x)f’(x)Sk(x)dx  + h l + (rg)j 


-xd/h 


< (C1(h,d)N(f(rg)',P)  + C2(h,d)N(frg*\P)]e 


(c)  Let  tfrg'/«H(x),  [frg$'/$]  (x)  and  [f'rg/$](x)  -*-0  as  x -*■  a 
and  as  x -*■  b along  T,  and  let  f(rg)",  f[2(rg)'$'  + rg$"]  and 
frg(*’)2  € B(P)  . Then 


(2.36) 


r(x)f"(x)Sk(x)dx 


- arg 

j—  J 1 


(rg)''  (Q)  [2(rg)^;*+  (rg)^'']  6kj 

Tj  5kj  + t]  h 


kJ  \| 
+ (r8)J^  ~J~  }| 


< (C1(h,d)N(f(rg)",P)  + C2(h,d)N(f (2(rg) '♦'+rg^"},P)  + C3(h,d)N(frg(4')/,P)] 


-xd/h 


Proof ; We  shall  only  prove  the  (b)-part  of  Theorem  2.9,  since  the  proofs 
of  the  (a)  and  (c)-parts  are  similar. 


We  find,  upon  integration  by  parts,  that 


(2.37) 


| r(x)f '(x)Sk(x)dx  - r(x)f (x)Sk(x) 

- | f (x)[r(x)s£(x)  + r' (x)Sk(x)]dx 


The  first  Cera  on  the  right -hand  side  vanishes,  by  assumption  of  the 
(b)-part  of  the  theorem,  while  by  expansion  of  the  second  part  of  (2.37), 
we  have 


(2.38) 


| r(x)f ' (x)S^(x)dx 


"l 


f (x) ( (rg) ' (x)S(k,h)e  $(x)  + (rg$') (x)S' (k,h).$(x) ]dx 


Hence  by  replacing  f In  (2.29)  by  the  integrand  on  the  right-hand  side 
of  (2.38),  and  noting  that  If  z £ dV,  then  |lm$(z)|  - d and 
|exp(^-  ♦(z)sgn  Im  $(z)]|  - e n<J/^  we  find  by  (2.29),  Lenina  2.9  and 
Theorem  2.3,  that 


I l 


+ h 


r(x)f '(x)Sk(x)dx 

r(rg)I 


< e 


[C.(h,d)|[f(rg)'](z)|  + C2(h,d)|(frg#')(«)|]|d«|, 
J aO  A L 


-nd/h 


which  Is  just  (2.35). 

Theorem  2.11:  Let  N be  a positive  Integer,  o a positive  constant, 
and  take  h ■ [nd/(aN)]^  . 

(a)  Under  the  assumptions  of  Theorem  2.9  (a). 


r(„)((„)Sk<„)«.,  I . 

*k  1 N* 


<2.39)  |fr 

where  depends  only  on  f,r,g,d  and  a; 


(b)  If  |lrgf](x)|  < exp[-a|*(x)|)  on  T,  Chen  under  the  assumptions 

of  Theorem  2.9  (b) , 

I f N f(rg)I  ) 1 

(2.40)  |jr  r(x)f(x)Sk(x)dx  + h ^ f,  {-jr*  + (rg)j  }| 


<.  K2  e 


-(irdaN) 


k - -N , -N+l N , 


where  depends  only  on  f,r,g,d  and  a; 

(•)  If  [f (2(rg) ' + rgi"/^' }] (x)  and  (rgf*'J(x)  are  bounded  by 
Kj  exp[-u|4(x) | 1 on  r,  then  under  the  assumptions  of  Theorem  2.9  (c) , 


<2-U)  \\t 


r(x)f"(x)Sk(x)dx 


- h 


«•  g<l>  fi(2) 

■?  f J™i  t«»  ; , (ra)  ^ || 

ji-K  “J  + ♦ 1-  h2  / 1 


< K3N%  e-("daN) 


k - -N.-N+l N , 


where  depends  only  on  f,r,g,d  and  a . 

Proof:  The  proof  Is  similar  to  that  of  Theorem  8.1  of  [13).  and  we  omit 
It. 

The  results  of  Theorem  2.10  are  especially  suited  to  the  solution  of 
linear  differencial  equations  via  a Galerkin  method,  for  which  the 
functions  (S^Jare  the  approximating  basis  functions.  We  remark  that 
we  could  have  obtained  alternate  expressions  of  j r(x)f ^ (x)Sk(x)dx, 

by  combining  equations  (2.29)  and  (2.27),  l.e..  If  rgf^  £ B(P),  then 
by  Eq.  (2.29) 


(2.42) 


If 


h rkgk  f(n)(iO 
r(x)f(n,(x)S.  (x)dx =- 

r k ’k 

< C1(h,d)M(rgf(n),P)e'1’d/h 


i 


and  we  could  now  use  (2.29)  to  approximate  f^(x^)  on  r . However, 
the  resulting  expressions  are  not  as  accurate  as  those  of  Theorem  2.10. 
Nevertheless  the  pair  of  equations  (2.27)  and  (2.29)  do  form  a powerful 
combination  for  purposes  of  solving  nonlinear  equations.  For  example, 
if<jG£B(P),  where  G ■ G(x,f  (x)  ,f ' (x))  , then 

G(xk,f(xk),f'(xk)) 


(2.43) 


H 


G(x,f  (x)  ,f  ’ (x))Sk(x)dx  - h 


♦'<V 


g<x.  ) 


< C1(h,d)N(G,0)  e 


-itd/h 


if  the  conditions  of  Theorem  2.7  are  satisfied  for  m - 1,  we  may  now 
replace  f'(xk)  *n  (2.43)  by  the  approximation 


N f 


;(1) 


(2.44) 


given  by  Eq.  (2.27) . 

The  approximating  expressions  of  Theorem  2.10  may  be  more  compactly 
expressed  by  means  of  matrices.  To  this  end,  let  m • 2N  + 1,  and  let 
S and  f be  column  vectors  defined  by 

•JB 


(2.45)  S_(x)  - 


■f-»  ■ 

S-N+1<X) 

f-N+l 

e 

• U ■ 

e 

e 

e 

e 

e 

V‘>  . 

.*»  . 

Corresponding  to  a function  u ■ u(x),  let  A (u)  denote  e dlegonel 


matrix,  whose  diagonal  elements  are  u(x_N>,  u(x_N+^) , . . . .uCx^)  and 
whose  off-diagonal  elements  are  zero.  Let  1^^  and  denote  the 

matrices 


(2.46)  I 


(1) 


0 -1 
1 0 

-1  > 

h.  _l_  — — 

2N  2N-1  2N-2 


1 

2 

-1 

0 


• • a 


1 

2N-3 


1_ 

2N 

1 

2N-1 

_1 

2N-2 


4s 


O) 

Jk 


(2.47)  I(2) 


JT 

3 


IT 

3 


(2N;  (2N-lj  (2N-2)' 


(2»>2 

2 

(2N-1): 


IT 

3 


4t 


With  this  notation,  EqJ>«.  (2.39),  (2.40)  and  (2.41)  take  the 
approximating  form 

| r(x)  f ' (x)S  (x)dx  - *,<'«)  U, 

* r 

f r(x)f"(x)S  (x)dx  - h[A  (^^)  ♦ £ lfl)  A (2(rg)' ♦ rg ♦"/♦') 


(2.48) 








__ 


Is  simply  related  to  a Toeplltz- 


I 


By  [4]  pp.  67-72,  the  matrix  I 


(1) 


m 


type  matrix,  by  considering  the  Fourier  series  expansion  of  -it  on  (-*,*). 

,<2) 


The  matrix  I 


is  a Toeplltz  matrix,  obtainable  by  considering  the 


Fourier  series  expansion  of  -t  on  [-s,ir]  . Thus  [4  p.65],  the  eigen- 
values of  are  iA^V 

m k 

while  the  eigenvalues  of  I 


k ■ -N,  -N+l , . . . , N , where  -it  < xj^  < it 

(2)  are  -X<2>,  k 
m k * 


0 < X(2)  < t,2 
k 


Indeed,  let  0 < X*2)  4 A*2^  < 


-N,  -N+1,...,N,  where 
...  < X<2>  < n2  . Then 

-n  — -«tj.  — “*  N 

2 

by  (4  pp.64  and  67],  since  x ^ 2 - 2 cos  x on  [ — ir , w ] , it  follows  that 
A(2)  >2-2  cos[ir/(2N+2)  ] - 4 sln2[  tt/(4N+4)  ] . That  is,  U£2)/X(2)|,  the 

— N N — N 

(2^  2 2 2 
condition  number  of  i''  is  bounded  by  x / { 4 sin  (ir/(4N+4)]}  (N+l) 

ID 


Summing  up,  we  have 
Theorem  2.12:  (a)  I 


(1) 


m 


k - -N,  -N+l N , (b)  I 


is  a skew-symmetric  matrix  having  determinant  zero. 

c 

(2) 


The  eigenvalues  ix£^  of  satisfy  the  Inequality  -it  < A^^  < i. 


values 


(£2\  where  4 sin2[ir/(4N+4) ] 4 A^2^  < it2 


is  a negative  definite  matrix  having  eigen- 

k - -N,  -N+l N . 


We  close  this  section  with  a derivation  of  the  formulas  of  Theorem 
2.10  for  the  case  of  the  important  intervals  [0,1],  [-1,1],  [0,»], 
and  [-«•,«]  . 


Ex.l: 


(2.49) 


[0,1] . In  this  case 


♦ (z)  - log(j^) 


♦ ’(z 


l V - (z:|arg(y^)|  < d}  . 

Let  us  assume  that  the  coefficients  r of  a second  order  equation  are 
analytic  and  bounded  in  V,  and  that  the  same  is  true  of  r'  and  r"  . 
It  is  then  convenient  to  take 


on  P,  and  if  on  10,1],  | f (x) | <,  C[x(l-x) ]°,  where  C end  a are 
poaitlve  constants.  If  f does  not  vanish  at  0 and  1,  we  replace 
f by  F in  the  differential  equation,  where 

(2.51)  F(x)  - f(x)  - a(l-x)  - bx 

and  where  a ■ f(0),  b ■ f(l)  . The  -^uncf/d* 


(2.52)  {sk(x))£B_N  - (x(l-x)S(k,h)oi(x)}k__N 


To 

x 


it  may  be  necessary  to  adjoin  1 - x 
is  unknown.  Differentiating  g and 


if  a is  unknown,  and 
♦ we  get 


(2.53) 


Hence 


(2.54) 


g'(x)  - 1 - 2x,  ,g"(x)  - -2,  ♦"(x)  - - 

x (1-x) 

(rg)(x)  - x(l-x)r(x)  (|$)(x)  ■ x2(l-x)2r(x) 

(-ffi -) (x)  ■ x(l-x)[x(l-x)r'(x)  + (l-2x)r(x) ] 

■ (-^ff^Kx)  - x(l-x)  [x(l-x)r"(x)  + 2(l-2x)r ' (x)  - 2r(x)] 
(X)  - 2x(l-x)r’(x)  + (l-2x)r(x) 
,(rgf')(x)  - r(x) 


Hence  we  gat  the  approximations  (2.48),  in  which  x^  ■ -j  + y tanh(kh/2)  . 
Ex. 2.  r ■ [-1,1].  In  this  case 





The  conditions  of  Theorems  2.9  and  2.10  are  satisfied  If  f Is  analytic 
and  bounded  on  P,  and  If  on  (-1,1),  |f(x)|  <,  C(l-x^)°,  where  C and 
a > 0 . If  f does  not  vanish  on  -1  and  1,  we  set  f ■ F + p In 
the  differential  equation,  where 

(2.57)  p(x)  . a l=x  + b i+* 


and  where  a - f(-l),  b ■ f ( 1)  , The  functions 

(2.58)  tSk(x))k-N  " {-|u-x2)S(k,hM(x)}“__N  ; 


to  it  u7  **e  necessary  to  adjoin  (l-x)/2  and/or  (l+x)/2  if  a 

and/or  b are  unknown.  Differentiating  g and  we  get 


f(rg)(x)  - ^(l-x2)r(x)  ; Oif)(x)  - (■^5-)2  r(x) 


(2.60) 


' 2 

2 


(-^^-Xx)  - <-^_>2r’(x)  - x(l=2L.)r(x) 


- (— 2r" (x)  - x(l-x2)r  ’ (x)  - (^)r(x) 


(2(rg)  ’ + rg^/VXx)  - (1-x  )r'  (x)  - xr(x) 
(rg4>')(x)  - r (x) 


Hence,  with  x^  - tanh(kh/2),  the  approximations  of  Theorem  2.10  take 
the  form  (2.48). 


Ex. 3:  The  case  T ■ [0,»]  . In  this  case 


(2.61)  +(z)  - log  z,  ♦ '(*)  - V ■ (z:  )arg  z\  < d)  . 


Suppose  that  the  coefficients  r and  the  derivatives  of  r are 


analytic  and  bounded  in  V . If  on  V,  |f(z)|  <,  C|z|°/(l+|z|) 


2a 


where  C,  a are  positive  then  it  is  convenient  to  take  g(z)  ■ z/(l+z)  , 


in  order  that  the  conditions  of  Theorems  2.9  and  2.10  are  satisfied. 

i®i/< . i i.2+a 


However,  if  |f(z)|  ^c|z|  /(1+|*|;  in  V,  where  C and  a are 
positive,  then  it  is  possible  to  choose  a simpler  form  for  g,  and 
Sk(x) , namely 


(2.62) 


g(x) 


Vm 


• X ; 


Sk(x)  - g(x)S(k,h)*4(x)  . 


In  this  latter  case 


FiJ.2.3 


•a. 


f(rg)(x)  - xr (x)  (|$)(x)  ■ xr(x) 

(^ffl  ) (x)  - x2r'(x)  + xr (x)  ; (x)  - x2r"(x)  + 2xr'(x) 

9 9 

(2.63) 

(2(rg) ' + rg$"/$')(x)  - 2xr'(x)  + r(x) 
l(rg4>')(x)  - r(x) 

kh 

The  approximations  now  take  the  form  (2.48),  in  which  x^  ■ e 

2 

If  f is  merely  bounded  on  V,  and  if  d ■ lim^ |m^x  f'(x),  then 
we  replace  f by  F in  the  differential  equation,  where 


(2.64) 


where 


(2.65) 


. , . , , a ,xb  . xc 

f(x)  - F(x)  + _ + _ + _ 


(1+x) 


a - f(0),  b - f(«),  c - b - a - d 


2 

If  the  limit  lim,  Nx  f*(x)  does  not  exist,  it  may  be  better  to  take 
vX***; 

2 

g(x)  ■ x/(l+x)  or  g(x)  ■ x/(l+x)  , depending  upon  the  problem. 


"T~ 


Ex. 4;  The  case  T ■ [-•,“]  . In  this  case 

(2.66)  ♦(*)-*,  ♦'(*)-!.  F,].2.+  Vofex.f 


and  V - P'  (see  Eq.  (2.12)).  If  the  coefficients  r of  the  differential 
a 

equation  are  analytic  and  bounded  in  , and  if  f £ B(P^) , we  simply 
take 

(2.67)  g'(x)  - 1,  (Sk(x))J-_H  - (S(k.h) (x) 

in  order  that  the  conditions  of  Theorem  2.9  become  satisfied,  and  provided 


that  f vanishes  at  ±°*  . The  conditions  of  Theorem  2.10  also  become 
satisfied  If  |f(x)|  ^ Ce  on  T . Then  “ kh,  and 

(rg)(x)  - r(x),  (^&)(x)  - r(x) , ) (x)  - r'(x)  ; 

9 9 

(2.68)  • (-^|^)(x)  - r"(x) , (2(rg) ' + rgf'/^H*)  - 2r'(x) 

(rg$')(x)  - r (x)  . 

The  approximating  equations  again  take  the  form  (2.48). 

If  f does  not  vanish  at  ±»,  we  replace  f by  F,  where 

(2.69)  F(x)  - f(x)  - _cx  (e"cxf(-«)  + eCxf(-)J 

e + e 

and  where  0 < c < ir/(2d)  . 

I 

' 

■ : 

I I 


3.  Examples  of  Applications 

In  this  section  we  shall  illustrate  the  application  of  the 
formulas  developed  in  Sec. 2,  on  the  solution  of  some  simple  ordinary 
and  partial  differential  equations. 

Ex. 1 : Consider  the  simple  problem 

(3.1)  fxx(x)  - -2,  0 < x < 1 ; f (0)  - f(l)  - 0 

This  has  the  solution  f(x)  ■ x(l-x)  . By  taking  r(x)  ■ 1 in  (2.54) 
and  combining  with  (2.48),  we  arrive  at  the  system  of  equations 

(3.2)  h[-2A  (x(l-x))  lf1}  A 6 - 2X;  + I^Jf, 

- -2hA  (x2(l-x)2)e 

T 

where  (1,1,. ..,1)  , T denoting  the  transpose.  Solving  this 

i U 11 

system  for  the  case  N ■ 16,  h -.7 SfH  ) x^  - tanh(kh/2),  we 

get  an  approximate  solution 
i6 

f(x)  * l f S(k,h)«<Kx)  (*(x)  - log[x/(l-x) ]) 

k*-16 

which  is  accurate  to  5 dec.  on  [0,1].  Similar  accuarcy  obtains  if  the 
-2  in  (3.1)  is  replaced  by  -2xa“2(l-x)0'2  , a,B>0  . 

Ex. 2.  f"  - f - f3/x3,  f(0)  - f(«)  - 0 . This  problem  was  solved 

kh 

by  different  procedures  in  [1]  and  [6].  By  taking  x^  • e and 
combining  (2.63)  and  (2.48),  we  get  the  approximating  system 

(3.3)  f I*1*  ♦ £ I<2)]f  - h A (x2)[f  - e ] 

L n +J&  — m '“-b 




n . ■*&*■■«***.  _ J 


-2  3 -2  3 -2  3 T 

where  e^  denotes  the  vector  [ x_Nf  x ^^f ••••»  *N  f^] 

The  solution  of  (3.3)  Involves  the  solution  of  a system  of  nonlinear 

j ^ j|| 

equations.  By  taking  h • .y/ty  >N  * 16  we  get  an  approximate  solution 
16 

f(x)  a l f.  S(k,h)o<t>(x)  ($(x)  - log  x) 

k— 16 

which  is  accurate  to  5 decimals  on  [0,®]  . 

Ex. 3.  u^  ■ u£  , 0 < x < 1 , t>0 

(3.4)  ^u(x,0)  » sin  xx,  u(0,t)  ■ u(l,t)  » 0 . 


In  order  to  get  zero  boundary  conditions,  we  set 
(3.5)  u ■ v + sin(irx)e  . 


This  yields  the  problem 

v^  - vt  - (n^-4)sin(irx)e-^t,  0 < x < 1,  t>0 

0.6)  f 

W(x,0)  - 0,  v(0,t)  - v (1 , t)  - 0 . 

We  solve  this  by  taking  our  approximating  basis  functions  to  be 

Sk(x)  - x(l-x)S(k,h)o6(x) , *(x)  - log[x/(l-x)) 

(3.7)  { * * * 

^S^(t)  - t S(£,s)»6  (t),  4 (t)  - log  t . 


The  problem  (3.6)  may  now  readily  be  reduced  to  a matrix  problem,  by 
proceeding  as  for  (3.2)  above.  Setting 


-x 


*The  solution  f satisfies  f(x)  ~ Ae  " as  x-?“  ,and  consequently.  It 

M A1 

may  be  necessary,  on  some  computers,  to  replace  by  ? where 


If  « 


in  order  to  avoid  underflow. 


‘H 


-A/ 


V 1 


1 


0 


(3.8) 


-N.-N+l 


V-N+l ,-N  V-N+l,-N+l  ...  v 


-N+l.Nl 


N.-N+l 


(3.9)  B - -2hA  (x(l-x);  + I^A  fl- X%)  + ± I(2)  (x  - i + i tanh  411)  •»  ') 

— m ' —in  *~m  h — m k 2 2 2 7 


1.1.  ,khv 


<3-10)  £ ■ ts£'m  '"/l  ] (Jn°)  " unit  matrix,  - ekS)  y 


(3.11)  D - hAm(x2(l-x)2)  , (xk  " \ + \ tanh(kh/2)  ) ; 

(3.12)  Z - sA^t2)  » (ck  “ eks  ) ; 


(3.13)  F - 


where 


-N,-N 


-N.-N+1 


f-n+i,-n  f-n+i,-n+i 


N.-N+l 


■ (x2-4)ain(irxk)e  * , 


-N+1,N 


we  arrive  at  the  matrix  eyatem 


(3.14) 


+ l/ce1  = F 


r 


Eq.  (3.14)  may  be  solved  by  diagonalizing  and  . If 

X_N,X_I|+1 . .Xjj  and  denote  the  eigenvalues  of  D 

and  CE  1 respectively,  obtained  by  taking  *~\p  &K  and  ZCt  'i  ' via 

e.g.  the  method  of  Golub  and  Relnsch  [3],  and  If  G - [g^]  “ x lpz  \ 

Y - [y^l  - X_1yz_1  then  y^  - 8^/(\  + W£) . «nd  V -JCYZ  . 

By  taking  h-.75/N^,  s-.S/N^  , N-16,  we  get  an  approximation 

16  * 
u(x,t)  * e sinnx  + £ v .S(k,h)o$(x)S(£,h)*$  (t) 

V,l— 16  U 

which  is  accurate  to  '1  dec.  on  [0,1]  * [0,®]  . 


Ex. 4 

(3.15) 


UXX  + Uyy  “ -1»  ^ S = l0'1!  X I0*1! 

u • 0 on  3S 


Letting  B and  D be  defined  as  In  (3.9)  and  (3.11)  we  now  get  the 


approximating  matrix  system 

(3.16)  C'J6Uf  U(D'Lb)r  - -W 

. a*  /♦*  • a 


where  U^-  [u^] , J1  - (h^J,  h^  - 1 


This  may  now  be  readily  solved 


via  the  dlagonallzatlon  of  pLQ  . By  taking  N ■ 16,  h *,75/H^ve  get  an 
approximate  solution 


16 

u(x,y)  * l u .S(k,h)o4(x)S(£,h)o4>(y) 

k,£— 16 


whlch  is  accurate  to  5 dec.  on  S . 


*•  ' 

4.  Error  Analysis 

For  sake  of  simplicity,  we  shall  restrict  ourselves  to  the  slaipler 
case  of  the  second  order  problem 

(4.1)  u"  + f(x,u)  - 0,  u(0)  » u(l)  ■ 0 

The  analysis  for  the  case  of  other  ordinary  or  partial  differential 
equations  is  somewhat  more  complicated,  but  may  be  carried  out  similarly. 

Throughout  this  section  a*ci »^2* * * * *^16  denote  positive  constants, 
and  h ■ (nd/(oN)^]  . 

In  the  notation  of  the  previous  sections,  we  take  $(z)  ■ log[z/(l-z) ] , 
and  we  take  tha  domain  of  analyticlty  to  be  V ■ (z:  |arg[z/{i-z)  ] | < d)  . 

! We  shall  assume  that  (4.1)  has  a (locally)  unique  solution  u^  which  is 

analytic  and  bounded  in  V and  which  satisfies  the  inequality 

1 

(4.2)  lu0(x)|  < CjX^l-*)8,  0<x<l, 


Definition  4.1.  Let  M(d,a)  denote  the  family  of  all  functions  v 
that  are  analytic  in  V,  such  that 


(4.3) 


v(0)  - v(l)  - 0 

gv"f  B(P),  |g(x)v"(x)|  < C2»a"1(l-x)°"1  on  (0,1); 

gf(*.v)£  B(P),  |g(x)f(x,v(x))  | < CjX^d-x)"'1  on  (0,1)  ; 


where 


(4.4) 


g(x)  - x(l-x)  . 


Ve  shall  also  assume  that  the  solution  of  the  Frechet  derivative 
problem 


*r  < • .-J 





(4.5) 


0"(x)  + fu(x,u(x))6(x)  - w(x),  6(0)  • 6(1)  ■ 0 


satisfies 


(4.6) 


|e(x)|  <.  c.| |a_1w| 


■ 


for  all  u£M(d,a)  such  that  | |u-Uq| | £ € where 
by 


(4.7) 

where 

(4.8) 


|f I | - sup  |f(x)| , 
xf(0,l) 


(A_1f)(x)  ■ -[  G(x,t)f(t)dt 


and  where  for  any  x & (0,1 J, 


(4.9) 


G(x,t) 


(l-x)t  If  0<  t <x 


x(l-t)  If  x<  t < 1 


Moreover,  we  shall  assume  that  If  ||u-u.||  <,  e,  then 


(4.10) 


||{A“1f(t,u(t))}||  < C5  . 


Is  defined 


Let  us  assume  that  we  have  found  an  approximate  solution 


N 


(4.11) 


u (x)  " I uw  S(k,h)e#(x)  (m  - 2N  + 1) 


k— N 


by  the  method  of  the  previous  sections,  and  let  us  set 


(4.12) 


6 ■ u - uA 

m m 0 


Then 

(4.13) 

e"(x)  + 
in 

fu(x,u(x) 

for  some 

u between  uA 

and  u , 

0 

m* 

(4.14) 

1 ».(«)! 

< C. 1 |u 
“ 4 ' ' m 

.-1. 


Now  by  Theorem  2.l|,  we  find,  by  taking  ^(x)  - g(x)S(k,h)e$(x) , 
x^  - j + j tanh(kh/2) , that 

f1  8<Xk) 

(4.15)  j tv"(x)  + f (x,v(x))]Sk(x)dx  » h lv"(*k)  + f^.v^))] 


and 


(4.16) 


fv" 

J0 


(x)Sfc(x)dx 


g"(x4) 


* h l v(Xj)[-,  6kj  + {2g,(Xj)+g(Xj)^"(Xj)/^,(Xj)} 


S(1) 

h 


d(2) 

+ g(x  )#'(x  ) 

3 3 n 


In  which  the  error  of  either  term  on  the  right-hand  side  of  (4.15)  is 

Jjt 

bounded  by  CgN~^e~^1,doiN^  and  the  error  of  the  right-hand  aide  of 

(4.16)  la  bounded  by  . By  our  process  of  solution,  the 

(1)  s(2) 

(4.17) 


numbers  uk  In  (4.11)  are  determined  such  that 


N g" 

h jL,  “l[*t  C + ,2,i  + ‘J  W ^ + «J*J 

•k 

♦ h ^ f(xfc,uk)  - 0,  k - -N,  -N+l N . 


Theorem  4.2:  Let  the  numbers  (k  ■ -N,-N+l N)  be  determined 

by  (A. 17),  and  let  ua(x)  be  defined  as  In  (4.11).  Then 

(4.18)  |ub(x)  - uQ(x)|  < C15N3/2e‘(l,daN)*,  0 < x < 1 . 


where  u^  Is  the  solution  of  (4.1). 


Proof : In  view  of  the  errors  in  the  approximations  (4.15)  and  (4.16), 
the  solution  of  (4.17)  is  equivalent  to  finding  a function  v6H(d,a), 


such  that 


(4.19) 


g(*k>  ek 

♦v(V)tv"(xk)  + f Uk,v(xk)]»  h~  ’ k * -N,-N+1,...,N, 


where  v(xk)  " u^»  an<*  where 


(4.20) 


lckl  < C8N*e_(,,dQ,N)  , k - -N,-N+1,...,N  . 


Since  v £ M(d,a),  it  follows,  for  any  t f (0,1),  that 


(4.21) 


(v"(t)  ♦ f(t,v(t))J  - l ^rj^-)Iv"(xk)  + f (xk,v(xk)))S(k,h)o4fr) 

. J y<;>CT.a|^ 


By  multiplying  (4.21)  by  4'(t)2,  taking  A 3 of  each  slda,  and  noting 
that  g(t)+'(t)  - 1,  we  gat 

v(x)  ♦ (A'lf(t,v(t)))(x)  - l £(v"  ♦ f(xk,v^jA'\#,(t)2S(k,h).6(t))(x) 

k— • fk 

f+ll)  - | 


Since  +'(t)  ■ 1 / [ t ( 1— t ) ] , It  follows,  by  taking 


t ■ [1  + tsnh(u/2)}/2,  x • [1  + tanh(w/2) )/2,  and  using  (4.8)  and 
(4.9),  that 


(4.23) 


Il(h,x)  = A 1{#’ (t)2ainc[(4(t)-kh}/h] )(x) 


-ft 


- tanh(w/2)  , ru-kh,, 

--tShg/a}  •lncf~h — ldu 


Th 


On  the  interval  [-•,«],  the  function  ll-tanh(w/2)]/(l-tanh(u/2)] 
Increases  monotonlcally  from  [l-tanh(w/2)]/2  to  1 while  on  tw,*»], 
the  function  (l+tanh(w/2)]/[l+tanh(u/2)]  decreases  monotonlcally 
from  1 to  [l+tanh(w/2)]/2  . For  this  reason,  it  may  be  shown  by 
a somewhat  lengthy,  but  simple  argument,  that 


(4.24) 


1 1.  (h,x)  | <,  4nh  . 


Similarly  if  x £ [0,1]  and  z £ 30,  we  can  show  that 

<*•«>  |X2(h.,)|  1 if 


since  Im  4(z)  - id  . 

By  means  of  (4.19),  (4.22)  and  (4.25),  Eq.  (4.21)  now  yields 


|v(x)  ♦ {A_r (t,v(t)))(x) | 

. . a kl  . 


(4.26)  4 j-K  ^ + Il1<h-*)l  |k|>H  ft*  * 

* ii2o>.«)i  j i -i 


Using  the  bounds  given  in  (4.20)  and  (4.24),  we  bound  the  first  sua 
on  the  right-hand  side  of  (4.26)  by  C9N3^2e~^1,dod^  ; using  (4.3)  and 

(4.24) ,  and  recalling  that  tanh(kh/2) , we  bound  the  second 

sum  on  the  right-hand  side  of  (4.26)  by  C^e  ' ' ; and  using 

(4.25)  and  the  fact  that  | sin[w$(i)/h] | ^ sinh[i»d/h]  if  z £ 3P, 
we  bound  the  Integral  term  on  the  right-hand  side  of  (4.24)  by 

2h  N(g[v"  + f(*  ,v)])/[dsinh(ird/h)]  ■ C^N  ^e  . Hence  for  all 

x £[0,1], 


(4.27)  |v(x)  + {A~1f (t,v(t)))(x) | < C12N3/2e“(ffdoN) 


Since  v 6M(d,a),  it  fellows  from  the  first  and  second  of  (4.3) 

that 

(4.28)  |v(x)  | < Cnxa(l-x)°,  0 < x < 1 . 

Furthermore,  since  v £-M(d,a),  and  since  ub  and  v coincide  at 

x_N,x_1)+3, . . . ,xN,  it  follows  that  [13,  Theorem  8.2]  for  all  x£  [0,1], 

(4.29)  |ub(x)  - v(x)|  < C14NV(ndaN)  . 

In  view  of  (4.5),  (4.6)  and  (4.10),  it  now  follows  that  for  all  x £[0,1], 

(4.30)  |u  (x)  - U_1  f(t,u  (t)))(x)| 

n in 

<,  |v(x)  + [A-1  f(t,v(t)))(x)|  C5|ub(x)  - v(x)| 

By  (4.14),  (4.27),  (4.29)  and  (4.30),  it  thus  follows  that  for  all 

* (■  10.1] 

(4.29)  |Bb(x)|  - |uB(x)  - uQ(x)|  < C15H3/2e-<1,daN)  . 


r 


This  completes  the  proof  of  Theorem  4.2. 

2 

Similarly,  It  may  be  shown  that  when  using  n ■ (2N+1)  points 


to  obtain  an  approximate  solution  of  a partial  differential  equation, 
such  as  (3.15),  the  error  Is  bounded  by  C.,!f^e  yN  < 5c  n^*e~Yn  . 

AO  — ^5 

2 

Indeed,  for  the  case  of  (3.15),  ve  may  take  C...  - 1 and  y ■ w . 
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